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Abstract 

The classical dynamics of a charged particle travelling in a laser field modelled by an elliptically 
polarized monochromatic electromagnetic plane wave is discussed within the time reparametriza- 
tion invariant form of the non-relativistic Hamilton-Jacobi theory. The exact parametric represen- 
tation for a particle's orbit in an arbitrary plane wave background beyond the dipole approximation 
and including effect of the magnetic field is derived. For an elliptically polarized monochromatic 
plane wave the particle's trajectory, as an explicit function of the laboratory frame's time, is given 
in terms of the Jacobian elliptic functions, whose modulus is proportional to the laser's intensity 
and depends on the polarization of radiation. It is shown that the system exposes the intensity 
duality, correspondence between the motion in the backgrounds with various intensities. In virtue 
of the modular properties of the Jacobian functions, by starting with the representative "funda- 
mental solution " and applying a certain modular transformations one can obtain the particle's 
orbit in the monochromatic plane wave background with arbitrarily prescribed characteristics. 

PACS numbers: 03.50.De, 52.20.Dq, 32.80.-t, 03.65. Sq 
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I. INTRODUCTION 



By the end of the 20th century Quantum Electrodynamics set an undoubted example 
of a theory allowing a thorough understanding and the ability to make predictions to an 
arbitrary precision the interaction of individual fundamental particles, photons and electrons. 
Until the development of the laser existing photon sources produced incoherent photons one 
at a time and the perturbative QED, based on an elementary single photon interacting 
with a single electron, was a perfectly adequate and accurate theory despite not taking 
into account the coherence aspect of the photons. Since the 1960's, after the advances in 
laser technology, the first theoretical studies of collective multi-photon relativistic dynamics 
began. Over the past fifteen years the revolutionary invention for a new technique of laser 
pulse amplification has led to the possibility to create the attosecond pulses of light beams 
of a very Mgh mten.ty and made the practical =tudy ot a new a.ea ot QED. ,Hys^cs „/ 
high-intensity laser-matter interactions possible (see e.g., textbooks 1] and recent surveys 

0.3). 

Contemporary renaissance of these research activities has opened up many new 
promising opportunities to develop our understanding of fundamental physics as well as the 
possibility of absolutely new applications. 

The necessity to modify the conventional single particle description of radiation-matter 
interaction at high intensities is arguably best illustrated by classical Thomson scattering, i.e. 
scattering of an electromagnetic wave by a free electron. The theory of ordinary Thomson 
scattering 0, 5] assumes that a linearly polarized monochromatic plane wave's phase ^ is 
approximated as uit — ki ■ x ^ uoit {dipole approximation) and the particle at rest responds 
only to the electric field 

E = Eq cosuit . 

Therefore the electron consequently executes simple harmonic motion at the same frequency 
as the wave, ul, along the electric field direction: 

g 

X = -EoCOSUJLt- (1) 



L 



The radiation from such an oscillating charge is calculated with the aid of the classical 



^ The monochromatic electromagnetic plane wave represents the simplest way to mathematically model a 
laser field, which is a reasonable assumption providing the transverse directions of the laser beams are 
much larger than the dimensions of the system considered. 
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Larmor radiation formula [5j which results in the universal, frequency independent, total 
cross section 



The assumption that the magnetic part of the Heaviside-Lorentz force of the plane wave, 

F = eE + -v X B, 

c 



can be ignored. 



V 

- I I-BqI I ^ I l-Eo 



c 

implies, taking into account the equality, ||-Eo|| = \\Bq\ \ , that the electron motion should 
be completely non-relativistic, f -C c. According to ([1]) the maximum electron velocity is 
Vmax = el IEqI l/mtUL and therefore the non-relativistic character of motion will hold true 



providing 



v'--=^^,E'o = -^h«l. (3) 



Here Xl is the radiation wavelength and the beam intensity II '■= cEq/Stt has been intro- 
duced. 

When the dimensionless intensity parameter rj is sufficiently large the dipole approxi- 
mation is no longer valid and magnetic force cannot justifiably be ignored. The electron 
motion becomes a nonlinear function of the driving force and relativistic effects will modify 
photon-electron scattering. The classical fully relativistic four-dimensional treatment of the 
light-electron scattering when the electromagnetic radiation is modelled by an electromag- 
netic monochromatic plane is well established. In this background the classical relativistic 
equation of motion for the electron with the full Heaviside-Lorentz force and taking into 
account the retardation effect can be solved exactly. This provides with an implicit para- 
metric representation for the electron trajectory^. In this solution the relativistic particle's 
position vector x is given as a function of the proper time which is in turn a function of x. 
The explicit solution x{t) is a subtle function of the physical time, known only in the form 
of an infinite series expansion over the harmonics with a fundamental frequency depending 



^ To the best of our knowledge J.Frenkel in 1925 was the fist one who presented the parametric relativistic 
solution Q . Since that publication this problem has been studied many times in different context (for the 
principal references see reviews and original papers Q) and became at issue for classical textbooks 

0, Si 01 ■ The brief history of strong-field physics written by the pioneer of this area can be found in the 
recent book 
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on the radiation intensity [12]. The appearance of the higher harmonics in the scattered 
radiation is a direct manifestation of this fact and results in the modification of the classical 
Thomson cross section ([2]) for a high-intensity laser beam ( see, for example, the basic ref- 



12 



13| and for the recent direct experimental observation of the second harmonics 



erences 

In order to fully understand the relativistic effects it is essential to have a clear under- 
standing of how the transition from non-relativistic motion to the relativistic motion occurs. 
Analysis (see for details 15|, ll6|) of the deviations from the non-relativistic regime shows the 
possibility to categorize the relativistic corrections into two groups; the first one is due to 
the influence of the magnetic part, v x B , of the Heaviside-Lorentz force and gives the first- 
order relativistic effects 0{v/c) . The second category are the "true relativistic" second-order 
effects 0(f^/c^) such as, for example, the well-known relativistic corrections to the kinetic 
energy, Darwin term and spin-orbit coupling. It was shown by H.Reiss 15|, ll6|, that there 
is an intermediate regime with a wide range of parameters (in terms of the laser's intensity 
& frequency) where the leading contribution to the deviation from the non-relativistic de- 
scription comes from the magnetic fields effects, while the "truly relativistic" f^/c^ effects 
can be still be neglected. This important observation gives justifies exploiting the combined 
methods, where the relativistic effects are taken into account only partially as corrections. 
Such an approach is, in particular, very attractive to the description of laser-atom interac- 
ions, where the whole formalism is intrinsically non-relativistic (see the recent publication 



17( 1 and references therein). It is also important since the widely used numerical technique 



becomes very cumbersome when passing to the fully four dimensional relativistic covariant 
description. 

Several partially relativistic approaches, taking into account magnetic field's influence on 



the particle's dynamics were developed (see, for example. [l5l. Il6l. Il8l| and references therein). 
The straightforward way to take into account the magnetic field is to go beyond the dipole 
approximation by including the full phase dependence uJit — ki-x in the corresponding plane 
wave potential^. Since the exact parametric solution to the relativistic Lorentz equation for 
a charged particle in plane wave background is known one could naively expect that the cor- 



^ Note however, the discussions in 19|,|20| of the possible artifacts arising from this way of partial relativistic 
consideration as well as necessity to include the radiation damping effects 



21 



22| 
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responding Newton equation can be solved exactly in an even more simple form. However, 
surprisingly, we were not able to find (in the extensive literature on laser-matter interactions) 
such an exact solution to Newton's equations of motion for a charged particle in plane wave 
background. In this note we aimed to fill this gap by providing the description of the inter- 
mediate region of the laser-particle interaction and to present an exact parametric solution 
to the non-relativistic equations of motion that is in a close analogy to the corresponding 
relativistic problem. As an application of the general formulae we will consider in detail 
the case of an elliptically polarized monochromatic plane wave. For this case we derive the 
explicit representation of a charged particle's orbit in terms of the laboratory frame's time. 
The solution is given in terms of the Jacobian elliptic functions, whose modulus depends 
on the background's radiation intensity and polarization. The particle's motion represents 
a drift displacement and an infinite sum of harmonic oscillations with the fundamental os- 
cillation frequency depending on the radiation intensity in a nonlinear way. Furthermore, 
owing to the modular properties of the Jacobian functions which is the existence of relations 
between elliptic functions with different periods, we attest an interesting duality between 
motions in backgrounds with a various intensities regimes. In particular, we show how a 
particle's trajectory in a monochromatic plane wave with a high intensity can be obtained 
from the "fundamental solution, " which describes the motion in a background with low 
intensity. 

The presentation of the material in the article is gathered as follows. In Section [III the 
non-relativistic motion of a charged particle in an external electromagnetic background is 
reformulated in a time reparametrization invariant fashion. The equal footing of time and 
space coordinates mimics the relativistic theory and enables us to use the conventional 
stationary Hamilton- Jacobi method to find a parametric solution to the non-relativistic 
equations of motion for a charged particle travelling in an arbitrary plane wave background. 
In Sections IIIII and IIVI we discuss in details a particle trajectories in the monochromatic 
plane wave background with an arbitrary elliptic polarization. In Appendix El we briefly 
comment on the instantaneous implementation of a gauge symmetry and Galilean boost 
transformation for a classical "non-relativistic" particle in an electromagnetic background. 
Appendix [HI sketches the derivation of a particle's trajectory for generic boundary conditions. 
Finally, in order to make article self-contained, Appendix[Clgives some mathematical features 
of the Jacobian elliptic functions used throughout the main text. 
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II. NON-RELATIVISTIC PARTICLE IN AN EXTERNAL FIELD 



A point non-relativistic particle, with mass m and electric charge — e, moving in an 
external electric field E and magnetic field B is infiuenced by the Heaviside-Lorentz force. 
A particle's trajectory x{t) may be determined from Newton's equations of motion 

m = e E{t, x{t)) + xB{t, a;(t)) . (4) 

The nonlinear equations (jlj) can be reproduced within the conventional variational principle 
of least action based on the following "non-relativistic" Lagrangian function'^ 

^ / da; \ m dx dx e dx ^ , , , , ^ , , , , 

if the external electric field E and magnetic field B are defined in terms of the gauge 
potential a?) = x) , a?)) in the standard way 

E{t, X) = x) - A{t, x) , (6) 

B{t,x) = Vx A{t,x) . (7) 

Here we intend to solve (jl]) for a special case of idealized laser field described by the so- 
called electromagnetic monochromatic plane wave. However, before restricting ourselves to 
this special case, we consider at first the more general plane wave background with a 
gauge potential of the form 

A,ix,t)=A,{0, (8) 
where is a 4-vector depending only on the light-cone coordinate 

n ■ X 

C = t , 9 

c 

with a unit constant 3-vector n pointing in the direction of a wave propagation. Also we 
assume that the Coulomb gauge is imposed which reduces to the condition 

n-A = Q. (10) 



The name "non-relativistic" is somewhat misleading because the Lagrangian (O is not Galilean invariant, 
one can speak only about an approximate Galilean symmetry for small particle velocities (see discussion 



23| and Appendix E|. 
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Bellow it will be shown that the Lagrangian system ([5]) with the plane wave background of 
type ([8]) is classically integrable and its solution can be represented in a parametric integral 
form. Furthermore, after specialising to the monochromatic plane wave background with 
an arbitrary elliptic polarization, we derive the explicit form of the trajectory in terms of 



the well-known Jacobian elliptic 
classical Hamilton- Jacobi theory 



unctions. To demonstrate this we exploit the ideas from 



24 



25 



26] as well as the Dirac constraint formalism \2l\ . 



A. The Dirac parametrization trick and the Hamilton-Jacobi equation 

Due to the explicit time dependence of the electromagnetic wave potential the Lagrangian 
([5]) describes a non- autonomous system. It is convenient to enlarge the corresponding con- 
figuration space in such a way that the resulting extended system becomes autonomous at 
the expense of being invariant under an arbitrary time parametrization. This method is 



often used in classical mechanics (see, for example, [25|], page 90 o r [26| . page 235) and also 
known to particle physicists as the Dirac "parametrization trick" 271. l28[|. 

The basic elements of this approach are the following. Starting from an arbitrary La- 
grangian system with Lagrangian C{x{t) ,da;/dt ,t) the configuration space is extended by 
considering time t as a new dynamical variable which, together with the other "spa- 
tial" coordinates x{s) , depends upon the auxiliary evolution parameter s . The dynamics of 
the extended system is determined from the degenerate, homogeneous (degree one), time- 
reparametrization invariant Lagrangian C* constructed from the initial C according to the 
rule: 

C (x{s) , t{s) , x{s) , i{s) \ := (^] C (xis) , ^/^ , tis)] . (11) 



\ds J \ ds ds 

Hereafter a "dot" over a variable denotes a derivative with respect to the evolution parameter 
s and we require that t{s) is a monotonic, increasing function of the new evolution parameter 
s 

dt , , 

Numerically the new classical action based on the extended Lagrangian (fTTl) is the same as 
the action of the initial system, however, it turns out to be invariant with respect to an 
arbitrary monotonic change of the evolution parameter 

s^s' = f{s). (13) 



8 



The extended phase of the system (fTT]) consists of 3 + 1 canonical pairs 



Zis) :-- 



t{s), Pt{s) 



P ■-- 



dC* 
dx{s) 



Pt 



dC* 
dt{s) 



(14) 



but due to the parametrization invariance f|T3|) the dynamics are constrained to develop on 
the surface of the phase space defined by the equation 



n:=pt + H, = 



(15) 



where He is the canonical Hamiltonian corresponding to the initial Lagrangian C . 

According to the Hamilton-Dirac description 27|, l28|] the constraint ( |T5l) plays a twofold 
role. First, it is a generator of the local symmetry transformation of the phase space coordi- 
nates f[T^ induced by the time reparametrization (1131) . Second, inasmuch as the canonical 
Hamiltonian derived from the Lagrangian C* is identically zero the dynamics are encoded 
in the constraint (ITS!) also. This constraint generates the evolution of the extended system 
via the Hamilton-Dirac equations 



z = x{s){z,n} 



(16) 



with an arbitrary function A(s) . The arbitrariness of the Lagrange multiplier A(s) reflects 
the freedom to use an arbitrary evolution parameter. Fixing it, by imposing an additional 
constraint 

Xis,x,t) = 0, 

allows one to show that any solution to f|T6l) either coincides with the classical trajectory 
of the initial Lagrangian C , if gauge y := t — s is chosen, or by using any other admissible 



gauge, is canonically equivalent to it 



28|. 



Applying this general scheme to the Lagrangian is transformed to 



(17) 



and the time reparametrization invariant Hamiltonian dynamics of the non-relativistic par- 
ticle is governed by the following Hamiltonian constraint 

2 



2m V c 



0. 



(18) 
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Our plan to find the trajectory of a non-relativistic particle in an electromagnetic plane 
background is the following. At first, exploiting the Hamilton- Jacobi method, we perform a 
canonical transformation to a free theory. Having the explicit form of this canonical trans- 
formation as well as the solution to the free equations of motion we will be able to write down 
the solution to the Hamilton-Dirac equation of motion (1161) . Then, fixing with a suitable 
gauge, the Lagrange multiplier function A(s) will be determined and finally the trajectories 
of a charged particle will be given in a parametric form analogous to the corresponding 
relativistic problem. 



B. Canonical transformation to a free system 



Following the basic idea of the Hamilton- Jacobi method consider a canonical transforma- 
tion 



x[s), p[s) 
t{s), Pt{s) 



Zo 



T{s), Ut{s) 



(19) 



that turns the constraint (1181) into the constraint of a free theory 



no = UT + 



2m 



0. 



(20) 



This canonical transformation "absorbs " the electromagnetic field and as a result the con- 
straint Tio generates according to ( |T6l) the simple free evolution: 



T(s) = To + / du X{u) , X{s) = Xo 
Jo 



— / du\(u) 



m 



(21) 



In this solution H-r and 11 denote constants of motion which contain all the information 
about the initial position and velocity of the particle when s = . The knowledge of the 
explicit form of the "absorbing transformation" Z ^ Zq is equivalent to solving the initial 
interacting problem and can be found in exceptional cases only. Fortuitously, it is the case 
for the problem we are considering here. 

The "absorbing" transformation (fT9|) can be established within the well-known method of 
generating function 2J, |25| using the S'2-function of old coordinates (t, x) and new momenta 
(Ht, n) of the form 

S2{t, x,UT,n) =tUT + x-n + ^(e, n) , (22) 
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with the unknown function 11) to be determined as follows. Write the old momenta pt 
and p as function of transformed coordinates 



Pt 



P 



dS2 
dt 



dS2 
dx 



^ nd^ 

c d^ 



(23) 



(24) 



Decompose all 3-vectors into their orthogonal and parallel components with respect to the 
direction of wave propagation, e.g., 11 = n_L + flu n, and by using the gauge condition 
ffTOj) we see that the constraint f|T8l) reduces to a free constraint fl20|) if the function JF is a 
solution to the equation 



1 d^ 

c de 



+ (mc — IIii 



(mc -n||)2 + vr(e,nx) 



(25) 



where 



W{i,Il^):=--Al + 2-A^-U^ + 2me^. (26) 

c 

The left hand side of ( !25|) is positive definite, therefore, a solution to (125!) is real function if 



(mc -n||)2 + vr(^,nx) > 0. 



(27) 



When the condition (ITTI) is satisfied and imposing the boundary condition J^{Q) = , we 
have the solution: 



J^(e,n) = -c(mc-n||)^ + c^ du sj{mc - Ii\\f + W{u, H^) . 



(28) 



Here we insist that the inequality (1271) is satisfied identically for all values of integration 
variable u. The importance of the restriction (1271) will be discussed later by analyzing the 
special case of a monochromatic wave background. 

Now substituting the function T from ( l28l) into ( l23l) and ( 12^ one can easily determine 
the expressions for a new momenta as function of the initial ones. The momentum which is 
canonically conjugated to the new time coordinate T is 



Pt-c 



mc — P\\) + \ {mc — p||)2 — , p±) 



and the new three dimensional momenta read 



n, 



p± , n,. 



mc + 



mc 



,)2-i^(e,px) 



(29) 



(30) 
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Using the generating equations 



T 



dS2 dJ^ 

^ = M = " + m 



(32) 

Il=Il{t,x,p) 



one can find the new coordinates as a function of the old ones. According to these results 
the time coordinate is unchanged 

T = t , (33) 
while the new three dimensional coordinates are 

e 

= x± + - r / duAJu), (34) 

X|| = X|| + + r[ du J inic ~p\\Y - Wiu, p±) . (35) 

\mc-p\i\ Jo V V 



Equations ( l30l) express three independent constants of motion, the two first constants 
coincide with the transversal momenta p± , while the third one Hn , can be interpreted as 
the longitudinal momenta of particle only in the asymptotic region where the interaction 
with the electromagnetic field is neglig ible^ Note also, from ([29]) and ([30]), it follows that 
the so-called "light-cone energy" represents a constant of motion 

— + P\\ = —— + nil = constant . (36) 
c c 

Suppose now that we are able to invert equations ( 1MI) -( 135I) . i.e. express the old canonical 
pairs {t,pt) and {x,p) as functions of (T,!!^) and (X , 11) . Further, since (T,!!^) and 
(X , n) are known from (!2T!) (up to a gauge fixing), these inverted expressions give the 
solution to the Hamilton-Dirac equations ([T6l) . To succeed in this inversion it is necessary 
to impose an appropriate gauge condition on the coordinates t and x. 

C. Light-cone gauge fixing and parametric solution 



The observation that the light cone energy ( 136]) is a constant of motion suggests a nat- 
ural gauge fixing condition. Namely, the evolution parameter s can be identified with the 



For further interpretation see equations ([51)1 and below. 
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canonical variable conjugated to the light cone energy 

X:=t{s)-^-s = 0. (37) 
c 

To find the Lagrange multiplier function A , note that according to the solution fl2Tl) 

T ^ = (1 h / duX(u), (38) 

c mc Jo 

where the initial conditions Tq and Xo have been set equal to zero for simplicity. Further- 
more, by using this relation and the equations fl33l) - fl35l) with the gauge fixing condition fl37]) 
implemented, we find that the Lagrange multiplier obeys the following integral relation 

mc f^"" 1 
duXiu) = — / du — , (39) 



With the aid of (139|) the equations (!33l) - (!35l) can be inverted with respect to the initial 
coordinates 

t{s) = mc du , , (40) 

h V(n|| -mc)2 + iy(M,n^) ^ ' 

r 1 

a^ll(s) = — cs + mc^ / d-u — , (41) 

" h V(n|| - mc)2 + iy(n, n^) ^ ^ 

,s nx--Ax(M) 

ic±(s) = c / dw , ^ ■ (42) 

Jo V(n|| - mc)2 + 1^(m, n^) ^ ^ 

The formulae f H0ll -( H2l) gives the parametric solution for a non-relativistic particle's trajec- 
tory in an arbitrary plane wave background. They are in a close analogy with the parametric 
solution of the corresponding relativistic problem [s], y, [lo| . 

In order to write down the three dimensional trajectory xii) , as a function of the physical 
time, it is necessary to invert fHOl) . i.e., to find s as function of t . It is not possible to 
write down an explicit formula for an arbitrary plane wave background. However, in the 
Section UTTl it will be shown how to solve this problem for a monochromatic plane wave. 



D. A particle orbit in a weak plane wave background 

Before considering the inversion of the integral (HOl) for a monochromatic plane wave we 
will sketch the possibility to solve the same problem for an arbitrarily "weak" plane wave 
in the form of and expansion in the intensity parameter. 

13 



According to (H0!) - (H2l) . if s as a function of t is known, s = f ^(t) , the classical trajectory 
can written in the form of the integral 

n I e ■ 



t- — [ dt'AM'\t')), (43) 
mc Jo 



m 
and 

zit) = ct - cf-\t) , (44) 

These formulae, with the function determined from (HOl) . gives a non-relativistic parti- 
cle's trajectory as function of LAB frame time in an arbitrary plane wave background. 

A "naive" non-relativistic limit of the solution for a particle's trajectory follows from (H3|) 
and (jUj) by assuming the validity of the formal 1/c expansion of the denominator of the 
integrand in the expression (HOl) . In this case, for small enough laser intensities, a closed 
form of the charged particle's classical trajectory, x{t), as function of the LAB frame time 
t can be written straightforwardly. 

Indeed, in the approximation mc — ~ mc and Ilj_/mc ~ , keeping only the leading 
term of the 1/mc-expansion of the denominator in ( l40l) we have 

t[s)=s + -rfl dua\{u). (45) 

2 .In 



In f H51) the normalized potential a± := (A^) with (...) denoting the time average 

and the dimensionless intensity parameter i], 

have been introduced. 

Therefore, for small intensities, the auxiliary time s in the leading rj order is 

1 /■* 

s = t rf I dua\{u) + ..., (47) 

2 Jo 

and the approximate form of a charged particle's trajectory reads 

r-t TT ft 



xAt) = 1 — cri I dua±_iu)^ rf / (lua?Au) + ..., (48) 

m Jo 2m Jo 



1 



z{t) = ^C7]' I duai{u) + ... . (49) 







The higher order corrections can be obtained in a similar way using, for example, the 



well-known Lagrange expansion method over the small parameter 
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29|. 



III. PARTICLE'S ORBIT IN A MONOCHROMATIC PLANE WAVE 

In this section we exploit the generic parametric representation for a particle's trajectory 
found above for the practically important case of a charge's propagation in the background 
of a monochromatic plane wave with an arbitrary polarization. 

We specify the gauge potential in equations fHOl) - (l42l) as 



Here, the four vector in brackets describes a monochromatic elliptically polarized electromag- 
netic plane wave with frequency ujl travelling in the 2;-direction. The parameter < e < 1 
measures the polarization in a way such that the boundary values e = and e = 1 cor- 
respond to linear polarization, while e = 1/v^ to circular polarization. In application to 
laser beams the profile function a{u) is usually assumed to be smooth and slowly varying 
(on the scale of oscillations) and vanishing at m ^ ±oo. In this article we are not intending 
to discuss a realistic laser and therefore, for the remainder of calculation, the pulse function 
is chosen to be constant a{u) := a. Formally this corresponds to a laser with an infinite 
length pulse. 

Any solution to the classical equation of motion for a particle depends on the laser field's 
characteristics as well as on the initial/boundary conditions on a particle's position and 
velocity. The laser field's characteristics used below are, the frequency ujl , the polarization 
e and the gauge invariant dimensionless intensity parameter ([3]) which for our choice of 
monochromatic wave potential is 



where now (...) denotes the time averaging over the period 2tt/uil- 

The dependence of a particle's orbits on the initial/boundary conditions is encoded via 
the first integrals 11^ and ny . The analysis of a generic boundary conditions is given in the 
Appendix [Bl It is shown there how one can invert equation (HOi) expressing the auxiliary 
evolution parameter s in terms of the physical time t via the Weierstrass elliptic function 
piujLt). Having this representation the remaining integrals fHT]) and P2|) with the arbitrary 
constants 11^ and Hy determine the explicit form of a particle's trajectory as function of the 
physical time t . However, to make the solution more transparent, we omit this rather tech- 




(50) 




(51) 
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nical work here and prefer to describe the orbits for a restricted but nevertheless informative 
initial conditions. 

Recall that the solution (Hni) - (H2I) is written when the initial conditions on the coordinates 
(in physical time t) read 

x{t = 0) = 0, (52) 

and the initial velocity v{0) := dx/dt{t = 0) is expressible via the dimensionless constants 
of motion 

/9+ := 1 - — , /3± = (A, P2) := — , (53) 
mc mc 

with the aid of relations 

v^{0) = c(3± -crie±, (54) 



vM = c-c^Pl-r]^e^ + 2r]e^-(3^ , (55) 

where for the choice fl50l) we have e± = {e, 0) . Since for small velocities the system possesses 
a Galilean symmetry, (see Appendix [A]) we can pass to a certain reference frame by specifying 
the constants v{0) . Below we restrict ourselves by insisting the vanishing of the transverse 
velocity^ 

(3± = Q. (56) 

Therefore, the orbits are specified by the constant, 

^ = ^ = 1-^(31-V'e\ (57) 

characterizing the particle's longitudinal velocity at t = . In order to have a real velocity 
we require > rj'^e^ . 

Now with this specific choice of constants, we shall find x{t), y{t) and z{t) as functions of 
the physical time t . Here we also note that due to the 2n- periodicity of the monochromatic 
gauge potential fl50|) all canonical coordinates will be treated as functions of a point on a 
circle with the standard trigonometric parametrization restricted to the principle domain 
[-7r/2,7r/2]. 



As it will be shown below the fixation f3± = , corresponds to zero average transverse velocity, {v±) = . 
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A. Fundamental domain and fundamental solution 

Now we want to describe a particle's trajectory as a function of the laboratory frame's time 
t. First of all let us preserve the condition fl7r|) which guaranties the monotonic character of 
the function t{s) (see equation (IT^ ). as well as reality of particle's trajectories is satisfied. 
For the monochromatic plane wave (!50|) the inequality (1271) with the vanishing transverse 
momentum, n_|_ = , can be rewritten as 

1 - /i^ sin^ u>0, (58) 

where 

We define three allowed domains as 

(I) < /i^ < 1 , (II) /i^ > 1 , (III) /i^ < . 

Below we show that having the knowledge of the solution for region (I), which we call the 
"fundamental domain", determines the solutions in all other regions. Namely, it will be 
demonstrated that all possible particle trajectories can be obtained from the "fundamental 
solutions" (solution depending on parameters from the "fundamental domain") by combi- 
nation of inversion /i ^ and rotation to the imaginary axis n t/j, . 
Besides this, we also present solutions for two special cases 

jj? = and fi^ = 1 . 

We emphasize here that these solutions can also be described by starting from the orbits 
from the fundamental domain and taking a corresponding limit. 

To prove the above statements we start with the analysis of the "fundamental solutions" 
and derive a particle's trajectory in terms of the physical time t . 
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1. Orbits for the fundamental domain, (I) : < /i^ < 1 . 



For the parameters from the "fundamental domain" (I) the inequahty (l58l) is true for all 
values of s from the interval — 7r/2 < n < 7r/2 . Equations fHUl) . fH2]) can be rewritten as 



t s] 



X[S\ 



y[-v 



UJLil-/3z) Jo 



du 



'1 — /i^ sin^ u 



ujl M 1- 2e' 



arcsm 



fi^ sm[ujLS) 



c 1 - £2 



ujl V l-2e' 



In 



/i^ cos(c<Jis) + a/1 — /i^ sin (cu^s 



1 + V/i' 



(60) 

(61) 
(62) 



and the component in the wave propagation direction is 



z(s) = ct(s) — cs . 



(63) 



This is the parametric solution to the equation of motion with parameter s from the principal 
interval 

-^<UJLS<-. (64) 

Now, thanks to L. Euler, A.M. Legendre, N.H. Abel and C.G.J. Jacobi, we know how to 
invert fl60l) . i.e. find the evolution parameter s as function of the physical time t. This can 
be done in terms of the well-known Jacobian amplitude function ( 29|, l30[ and Appendix [C]) 



uls = am(c<j^t, /i) 



with modulus := y ii^ and the non-relativistically Doppler shifted frequency 



u'l := iOL (1 - Pz) ■ 



(65) 



(66) 



For the values of s from the interval (1641) the amplitude function is a well defined increasing 
function defined on the interval 



- K(/.) < uj'^t < K(/i) 



(67) 



where K is the "real" quarter period of the Jacobian elliptic functions, (1C7I) . Therefore, one 
can consider the transformation from the evolution parameter s to time t as well-defined 
change of coordinates on a circle. 
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Substituting the expression for the evolution parameter in terms of the physical time (165|) 
into fl6T]) . fl62l) and fl63l) and using the properties of the Jacobian functions we arrive at the 
representation of the classical trajectory 



xpit) = 



1 -2e^ 



arcsin [/isn(co'^t, /i)] 



c 



In 



/icn(c<j^t, /i) + dn^cu^t, /i) 

1 +/U 



and 



Zpit) = ct 



am 



(68) 
(69) 

(70) 



Here the subscript "F" is written to emphasize that ( I68l) -( I70|) corresponds to the trajectories 
with the modulus from the fundamental interval < < 1 . 

Now we shall consider all the other domains of parameters and, using the properties of 
the Jacobian functions, will give explicit expressions for their corresponding trajectories. 



2. Orbits for the second domain, (II) : fj? > 1 . 

When analyzing this domain two peculiarities should be taken into account. First of all, 
when /i^ > 1 the inequality (158|1 is true only if 

sin^ M < /i^ , /i := — . (71) 
- 

This means that the upper limit to the integral in fl^Ul) lies in the interval 

— < — arcsin (/x) < ujls < arcsin (/i) < ^ . (72) 

Second, since the standard integral representation of the amplitude function ( ICip is 
defined for a modulus from the fundamental interval, < /i^ < 1 , some simple mathematical 
manipulations are required to rewrite the integral (!60|) in such a form that its integrand 
depends on the inverse modulus, fi , instead of modulus fi . Namely one can easily verify 
that pUj) can be written as 

-I /•arcsin(^ sin(a;j^s)) -i 

tis) = ^ du . (73) 

Therefore, the relationship between the evolution parameter s and time t is 

ujls = arcsin (//sn (u;^/xt, fi)) . (74) 
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When the parameter s is contained in the interval (172|) relation (17^ defines an increasing 
function on the interval 

-K(/i) <cj>t <K(/i). (75) 
Finally using (ITill and relations ( 1C2I) we have 



x(t) 



V 1 - 2 £2 



1-^2 



In 



am(u;^/it, /i) 



cn(co'2/it, /i) + /i dn(ci;2/it, yu) 
/i+ 1 



and 



z(t) = ct arcsin (/isn(ci;^/it, /i)) . 



(76) 
(77) 

(78) 



These formulae describe a particle's trajectory when the parameter /j,^ takes its value in the 
second domain, /x^ > 1 . 



3. Orbits for the third domain, (III) : /i^ < . 



If /i^ < the inequality ( l58ll is true for the whole principal interval [— 7r/2,7r/2] . Now 
introducing the positive definite parameter > ,/i2 := — , the relation (l60|l reads 



t(.) 



1 



dw 



while the solution for the spatial coordinates is 

c rje 



(79) 



x{s) 



arsinh [k sin(ci;i^s)] 



In 



IK 



K cos{ujls) + z V 1 + ^2 sin {ujls) 
i + K 



(80) 
(81) 



Now again a change of integration variable must be done in (179|) in order to get the integral 
expressible in terms of the Jacobian amplitude with modulus from the open interval (0, 1), 
thereby guaranteing that all functions are single- valued and continuous. It is straightforward 
to check that ( !79l) is equivalent to 



t{s) 



0(s) 



du 



1 K sin u 

k'2 



(82) 
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where the upper hmit of the integral is 



K Sim LOS 



arcsm 



a/1 + sin^(co's) 



(83) 



We have now achieved the goal that the modulus of k/k' G (0, 1) , for all values yU^ < 
Therefore, the inverse to fl82l) reads 



0(s) = ajuiuj'j^Kt, k/k') , (84) 
or directly for the evolution parameter we get 

With the aid of (!85|) one can finally rewrite the particle's trajectory in terms of the Jacobian 
elliptic function with modulus k/k' 



/ \ c rie . , 
x{t) = arsmh 



y{t) 

z{t) 



C rjy/T 



K sn {uj'j^n't, n/ tx) 
K,' dn {uj'^n't, k/k,') 

1 — m cn {uj'j^k' t, k/ k') 



In 



IK 



ct arcsin 



1 — m)dn {uj'i^K't, k/ k') 
1 sn (cu^K't, k/ k') 



(86) 
(87) 
(88) 



k' dn {uj'^K't, k/ k') 

Equations fl86|) - fl88|) describe a particle's trajectory in the background characterized by 
/i^ < , and are well defined on the interval 



K{k/k) < io'^n't < K{k/k) . 



(89) 



4- Degenerate orbits, = &: /i^ = 1 . 

Here we consider the two special cases remaining from the previous considerations. 
Namely, the special orbits with the parameter /x^ = and /i^ = 1 will be presented. Note 
that the vanishing modulus corresponds to a special circularly polarized monochromatic 
plane wave, = 1/2 or the trivial zero background case. 

If /i^ = 0, then the parametric solution PUI) determining the time t as a function of an 
auxiliary parameter s takes the simple form 

t(s) = ^^s, (90) 
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while the equations (H2!) for the spatial components orthogonal to the direction of the wave 
propagation reduce to 

1 C , . I- C 2 {^LS\ , ^ 

x{s) = -= — r/sina;LS, ?/(s) = - V2 — sm ( — l- (91) 

Therefore, it follows from (Mil) , that for the degenerate case of zero modulus /z^ = the 
trajectory is 



x[t) = —ii—smujrt, vit) = =r7— sm 



z{t) = PzCt. 



From these expressions we see that for a circular polarized monochromatic wave all nonlinear 
effects disappear and one can choose such a reference frame where the particle's motion 
appears as a pure harmonic with the non-relativistically Doppler shifted laser frequency cj^ . 
For the orbit characterized by /i^ = 1, the equation fHOl) gives 

t{s) = — artanh (sinw^s) , (92) 



and from fl42l) it follows that 



^is) = -\l ^_2^2 = Y 1 _ 2 g2 ~ ^'^ (^"^ ^ls) ■ (93) 

This finalizes our consideration of all possible particle's trajectories in a generic monochro- 
matic plane wave. 

In the next subsection we will briefly outline how all possible particle's trajectories can 
be categorized into the equivalent classes with respect to the action of the SL{2, Z)/(l, —1) 
group. 



B. Modular properties of orbits 

As shown above, the nonlinear dependence on the polarization and the intensity of the 
radiation background of a particle's trajectory is encoded in the modulus of the elliptic Jaco- 
bian functions. The doubly periodic elliptic functions have a remarkable property related to 
a certain symmetry of modulus transformations. Namely, the elliptic function with periods 
Wi and W2 can be algebraically expressed through another elliptic functions with periods w[ 
and W2 if periods are related by the so-called unimodular transformations 

w'l \ I a b \ I wi 
w'2 j \cdl\w2 
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where the entries of the 2x2 matrix are integers a,b,c,d G Z and ad — be = 1 . The 
transformations fl94|) may also be treated also as a subgroup of the Mobius transformations 
in the upper half of the complex r plane into itself 

T = —-7-, r := — . 95) 

a + OT Wi 

For the problem we are dealing with here this modular equivalence exposes in the intensity 
duality, i.e. there exist a specific correspondence between the motion in backgrounds with 
various intensity regimes. Any trajectory with an arbitrarily prescribed intensity can be 
connected to the solution from the fundamental domain with the aid of a certain modular 
transformation. Particularly, if we assign to a particle's trajectory with low intensity a 
certain modular parameter r, then the trajectory for the high intensity conditions are related 
by r ^ r/(l + r) . Indeed, using the relations between the Jacobian functions whose moduli 
are inverse to each other, (cf. formulae flC27|) ). one can verify that the expressions (I76l) - fl78|) 



with /i^ > 1 follow from the fundamental solution (l68l) - (|70l) by the substitution /i ^ 1//^ • 

x{t I /i) = xp{t \-), (96) 

or in terms of the modular parameter r 

x{t\r) = XF{t\-^). (97) 
1 + r 

Analogously, if /i^ < 0, the trajectories (!86l) - (!88l) are connected to the fundamental solution 
by the shift transformation flC26l) . 

x{t\T) = XF{t\l + T) . (98) 

Note also that the special trajectories with //^ = and //^ = 1 , considered above, coincides 
with the corresponding limits of the fundamental solution taking into account that the 
Jacobian functions are degenerate to the trigonometric (]C13I) and hyperbolic functions (1C12P 
for moduli /i = and /i = 1 , respectively. 



IV. ANALYSIS OF THE PARTICLE'S TRAJECTORY 



Now we analyze in greater detail the trajectory with parameters from the fundamental 
domain and clarify some of the physical features of this solution. 
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From the solution (l68!) - (170|) and (1C4P - (1C6P one can derive an expression for the particle's 
velocity 



Vy{t) 

Vz{t) 



-cqe cn (cj^t, /i) 



— cr^Vl — £^ sn {uj'it, fi) , 
c- c{l- P^)dn{uj'^t, n) 



(99) 

(100) 
(101) 



Periodic properties of the Jacobian function (see eqs. (ICQO in the Appendix [C]) tell us 
that the components of the charged particle's velocity in the plane orthogonal to the wave 
propagation are periodic functions of time with period 



4K 



UJ'r 



271 
Up 



(102) 



while in the direction of propagation the oscillation's period is half the size Tp/2. The 
fundamental circular frequency of the particle's motion, up, differs from the frequency of 
the laser field 



UJp 



2K 



(103) 



From this expression we see that, apart from the pure kinematical non-relativistic Doppler 
shift (!66|) . the frequency of a particle's oscillation depends on the laser's intensity through 
the real quarter period K flC7|) of the Jacobian functions. The presence of K in fllOSp 
exposes a new property of a particle's dynamics which is beyond the dipole approximation. 
A particle oscillates at frequency that depends on the laser intensity and polarization in a 
nonlinear way. For the low intensity regime, -C 1 , the period of a particle's oscillation can 
be represent in the form of an expansion (ICISP as 



27r 



1 + 



l-2e^ 



■T] + 



1-3 

2^4 



:i-26 

(1-/5. 



2\2 



(104) 



It is well-known that, in contrast to the dipole approximation where the particle's motion 



is a pure harmonic, the relativistic dynamics ex 
the particle which depends on the intensity 12, 



libits a drift action of the laser field on 
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32]. In our consideration, when the 



relativistic effects are partially taken into account, this effect also can be seen. For the 
non-relativistic case with our choice of the initial conditions, 11^ = , the mean velocity for 
the transverse direction vanishes 

{vj_) = 0. (105) 
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While the drift in the direction of propagation is a nonhnear function of the laser beam's 
intensity 

{v.) = c-^{l-f3,). (106) 
This drift velocity (vz) for small intensities at leading order reads 

(v^) = vM (l - 4^rz^ V'^+---- (107) 

Another new feature when comparing to the the dipole-approximation is the appearance 
of the higher harmonics in the particle's motion. This can be seen from our solution with 
the aid of the well-known Fourier series expansion for the Jacobian function [29|, 130|. Using 
these formulae collected in Appendix [Cl one can represent the trajectory as 

Ace ^ o"-^/2 
x(t) = > \- T^—-sm(2n-l)ujpt, (108) 



yit) = — V- ^- — — sinMn-- u;pt, (109) 



z{t) = {v,)t-—f2 /^f sm2nijpt, (llO) 

, ^(1 + q ) 

n=l ^ ' 

where q is the so-called nome parameter 

q := exp (-'^^j ■ (HI) 



Note (see eq. flC22p ) that the nome q for small intensities is approximately 



When the intensity parameter is small one can perform a Galilean boost with the velocity 
V := — (0, 0, (vz)) to the so-called average rest frame (ARF), frame where the mean parti- 
cle's velocity vanishes. In this frame the particle's motion represents only the superposition 
of harmonic oscillations with the fundamental frequency up . 

In the ARF frame equations fllOSp - flllOl) reduce, for small intensities, to the following 
expressions in the leading order of the r^-expansion 

CB 

a^ARpW = r^sina;^t, (113) 

VA^At) = - ^^^^HZ ^ (1 _ cos^^t) , (114) 
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Pictorially these formulae describe the orbits shaped hke a figure of eight. This form of the 
orbit is well-known from the parametric relativistic solution j8l.l9l.flol and (11131) - flllSp are the 
small intensity approximation to the explicit representation of the relativistic trajectories. 
We omit here the detailed comparison with the relativistic solution since it requires the 
knowledge of a relativistic particle's trajectory as a function of the physical time as well as 
a careful analysis of the frame dependence of our partially relativistic solution. We intend 
to do this in a forthcoming publication. 

All mentioned features, the Doppler shift, dependence of the particle's oscillation fre- 
quency on the laser beam's intensity as well as the presence of higher harmonics in the 
particle's motion, lead to a several important phenomena. Among them there are the non- 
linear modification to the classical Thomson scattering and a charged particle's mass/energy 
shift in the electromagnetic background radiation. 



V. CONCLUDING REMARKS 



Most textbooks on the classical dynamics of particles pay tribute to a historical tradition. 
They begin with the consideration of non-relativistic mechanics and, after introducing the 
basics of relativity, discuss the specific features of the corresponding relativistic problem. 
The exception to this rule is the charged particle's dynamics in the background of an elec- 
tromagnetic plane wave. All well-known sources either discuss the problem in the dipole 
approximation {4] or solve the problem in the framework of relativistic mechanics directly 



io|. 



In the present note we "restore historical justice" by considering the classical problem 
of a Newtonian particle's motion in a given electromagnetic plane wave field beyond the 
above mentioned dipole approximation. We obtain an exact representation for a particle's 
orbit in the parametric form which is analogous to the well-known relativistic solution. 
Furthermore, we show that the three-dimensional non-relativistic orbit of a particle's motion 
in a monochromatic arbitrarily polarized plane wave admits, as an explicit function of the 
laboratory frame's time, the exact solution in terms of the doubly periodic elliptic functions. 
This is in contrast to the relativistic trajectory which is known explicitly only in a four 
dimensional parametric form. The derived solution explicitly exposes the presence of higher 
fundamental harmonics in a charged particle's motion as well as nonlinear dependence of the 



26 



oscillations on the intensity and polarization of the monochromatic background. Another 
interesting characteristic of the trajectory we discovered is an intensity duality between 
motion in different radiation backgrounds. Since a given intensity defines the modulus of 
the Jacobian functions, a particle's trajectory in backgrounds whose intensities are connected 
by the modular transformations are simply related. 

Finally it is worth commenting on the sinuous way to derive the exact particle trajectories 
within the reparametrization invariant form of the Hamilton- Jacobi method presented here. 
We prefer to follow it since this approach incorporates both a classical as well as a quantum 
treatment. The knowledge of the classical generating function connecting the interacting and 
free systems helps to construct the corresponding quasiclassical description and in particular 
can be exploited for study of the laser-atom interactions Furthermore, the suggested 
approach mimics the four dimensional relativistic treatment which is useful to study the 
deviation from non-relativistic motion. All these issues we intend to discuss in detail in future 
publications. Apart from this, based on the explicit form of the non-relativistic trajectory 
we plan to study several important effects, including the non-linear Thomson/Compton 
scattering, electromagnetic dressing, charged particle acceleration by a laser field etc., in the 
transition regimes between the non-relativistic and relativistic cases. 
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APPENDIX A: REMARKS ON "NON-RELATIVISTIC" SYMMETRIES 

Below we discuss the realization of the gauge symmetry as well as the implementation of 
the Galilean boost for a charged "non-relativistic" particle travelling in an electromagnetic 
background. 
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1. Gauge transformations 



In terms of the 3-dimensional notations for a gauge potential = , A) and coordi- 
nates = {ct , X, y, z) = {ct , x), the gauge transformation reads 

d 

A{t,x) A'{t,x) = A{t,x) + —n(t,x) , (Al) 

ox 

x) ^ x) = <l>{t, x)--^^^ n{t, x) . (A2) 

Under the transformations (lAip and (]A2p the "non-relativistic" Lagrangian (jS]) is not in- 
variant: 

^ ( Ax \ ^ ( Ax \ ed^, ,,, 

Here d/dt denotes the total derivative 

A Ax d d . 
di d^ ■ 9t ■ ^ ^ 

However, the variation flA3l) being the total derivative of the function e/cfi(t , a3(t)) , does 
not affect to the classical equations of motion. 

The changes of gauge potentials ( lAll) and (]A2[) are canonical transformations with the 
generating function 

F{x,P,t) = x-P -n{t,x{t)). (A5) 

2. Galilean boosts 



Under a Lorentz boost in the x direction with the factor 7 := 1/a/1 — jc? 

t' = 7 - ^ , z' = 7 (x - , = 1/ , = z , (A6) 

a gauge potential = ($, A) being a Lorentz 4- vector transforms as 

<f' = 7($-^^x) , A/ = 7(^x-^'^') , V = ^?;, AJ = A,. (A7) 

For small velocities, neglecting terms of order V^/c^ and higher the Lorentz boost (IA6p 
reduces to the Galilean boost 

t' = t , x =x-Vt, y' = y, z' = z, (A8) 
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while the gauge potential changes as 



= $ A. 



■X ) 



A. 



■X 




— ■ 



(A9) 



c 



c 



With the help of these relations, keeping only the leading terms in V/c, we have 



C^C = C- (vx-lvh 




(AlO) 



note that V is not infinitesimally small. The Lagrangian is invariant under the Galilean 
boost up to a total derivative only. 

Both variations of the Lagrangian, under the gauge transformations (1A3|) as well as under 
the Galilean boost flA10|) . are not important at the classical level but can lead in general 
to a nontrivial quantum phenomena. Particularly, they lead to the appearance of the so- 
called 1-cocycle for the wave function. The action of the unitary operator Ug generating the 
Galilean boost on the wave function reads 



It is interesting that for the charged particle this cocycle can be trivialized. Indeed, if the 
Galilean boost is accompanied by the gauge transformation generated by flA5|) with the 
special gauge function 



then the Lagrangian ([5]) remains unchanged and cocycle is removed. 
APPENDIX B: GENERIC BOUNDARY CONDITIONS 

Here we briefly state how one can express the physical time t in terms of the auxiliary 
evolution parameter when the generic boundary conditions on a particle's position and 
velocity are imposed. To find this dependance one can proceed as follows. Plugging the 
expression for the gauge potential fl50|) into the function W{u, from fl26|) we have 

W{u, Il±) = —1] m c [e cos m + (1 — e ) sin u) + 2r]mc (Hi cos m + 112 sin u) . 

The integral in (HO!) with this expression represents the so-called elliptical integral. Indeed, 
using the universal trigonometric substitution x = tan(u/2) it can be rewritten in the 




(All) 
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Weierstrass form: 



2 r-tim{u)Ls/2) 

t{s) = — / dx — -== 

ujl Jo y^fix) 



with the fourth order polynomial 
whose coefficients are 



(Bl) 



(B2) 



(B3) 
(B4) 



According to the classical result (see e.g. 0]) attributed to Weierstrass, the integral (IBT 
can be inverted: 



tan j 



V/(o).'(fO-5 






+ ^/(0)/"'(0) 


2 




4/"(0) 





(B5) 



Here the number of primes over the polynomial ( ]B2p denotes the order of the derivatives 
with respect to x. In (IB5I) the Weierstrass doubly periodic function p{z] g2 , gz) depends on 
two invariants g2 and (73 of the polynomial (]B2p 



^2 : 
^3 : 
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These equations show that the dependence of any particle's trajectory on the initial con- 
ditions is accumulated in the invariants 5'2,5'3 and functionally is quite subtle. Particulary 
there is no rotational degeneracy with respect to the vector n_L . 



APPENDIX C: VOCABULARY ON THE JACOBIAN ELLIPTIC FUNCTIONS 



Following the classical textbooks 
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30| and handbook [33| we collect the basic formulae 



on Jacobian elliptic functions which are extensively used in the main text. 
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The amplitude function := am(z,/i) is the inverse of the function defined by the 
integral 

^(0) = r d^^=^—= . (ci) 

Jo v 1 — /i^ sm^ 

The amphtude ain{z , /i) is an infinitely many-valued function whose principal domain of 
definition for real z is (— K , K) . The value of this constant K is determined by the modulus 
fi via the so-called complete elliptic integral, see ( ]C7p below. 

Three basic Jacobian elliptic functions, sn{z,fi), cn{z,fi) and dn{z,fi) are analytical func- 
tions of the complex variable z everywhere except at the simple poles. These elliptic function 
are expressible in terms of the amplitude function 

sn(z, /i) = sin (am(z , /i)) , cn{z, fi) = cos [am{z , fi)) , — am{z, jj) = dn{z, fx) , (C2) 

and satisfy the basic algebraic 

sn^2; + cn^2; = 1 , /x^sn^^; -|- dn^2; = 1 , (C3) 

and differential relations 

-^sn(z, /i) = cn{z , n) dii{z , n) , (C4) 
az 

d 

— cn(z, yu) = —sn(z,fi)dn(z,fi), (C5) 
dz 

— dn(z, /i) = — /i^sn(2; , /i) cn(z , yu) , (C6) 

which show the analogy of the Jacobian functions to the trigonometric functions. 

Functions sn{z,fi), cia{z,fi) and dia{z,fi) are doubly periodic functions of z. Periods of 
sn(^, fi) are 4K and 2iK' while periods of cn{z, /i) are 4K and 2K + 2iK' . Function dn(2;, fi) 
has periods 2K and 4zK' . The "real" (K) and "imaginary" (K') quarter periods are real 
numbers given by the complete elliptic integrals 

K(^) := \ , , (C7) 

Jo V 1 — /i^ sm V 

'•■k/2 



^K'(/z) := ^ f d^ ^ =. (a 

Jo J I - il - u?) sviT d 



1 - (1 -/i2)sin'^9 

The Jacobian functions as functions of the modulus are single valued on the complex plane 
with two cuts [1, oo) and (— oo, 0] . 
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The discrete symmetry. The Jacobian functions cn and dn are even functions , while sn 
is odd. They obey the relations 

sn(u + 2mK + 2niK' = sn(ti , //) , (C9) 

cn(w + 2mK + 2mK' , /i) = (-)™+" cn(M , /x) , (CIO) 
dn(M + 2mK + 2mK',/x) = dn(M , /x) , (Cll) 

where n,m E . 

Two degenerate cases. The doubly periodic Jacobian functions degenerate to other func- 
tions when one of the periods become infinite, that if is or 1 . 

When II — 1 the real quarter-period K = oo and Jacobian function degenerate to the 
hyperbolic functions 

1 1 
sn(u, 1) = tanhu , cn(u, 1) = — - — , dn(u, 1) = — - — . (C12) 

coshw coshu 

li II — the real quarter-period is finite, K = 7r/2 , but the imaginary quarter-period is 
infinite and the Jacobian function degenerate to the trigonometric functions 

sn(M, 0) = sinii , cn{u,0) — cosu , dn(M, 0) = 1 , (C13) 

Small modulus expansions. When modulus is small enough the Jacobian functions can be 
approximated as 

am(2;, ii) = z — -ii'^ [z — sin(2;) cos(^)] -|- o(/^^) , (C14) 
sn(2;, n) — sin{z) — ^n'^[z — sm{z) cos(2;)] cos(2;) -|- o{n^) , (C15) 
cn(^, n) = cos(2;) + ^11^ [z — sm{z) cos(2;)] cos(2;) -|- o{n'^) , (C16) 
dn{z, 11) = 1 - ^ //^ sm^{z) + o(/i^) . (C17) 

The quarter period K has a small modulus expansion of the form 



(C18) 
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The Fourier series expansions: These are 



smz. 



cmz 



2vr ^ g"-V2 

^ sin(2n-l)— , 



/iK ^ I- q 

71=1 

27r g 



71=1 

dn(.,;.) = ^ + ^LrT?^-^2n- 

n=l 



(C19) 
(C20) 
(C21) 



with the so-called nome, or Jacobi parameter q = exp (— ttK'/K) . Similarly to the quarter 
periods the nome q can be expanded in powers of the modulus 



(C22) 



The modular transformations. The ratio r := iK' /K serves as an important modular 
parameter''. Under the unimodular transformation 

c + dr 



a,b, c, d G Tj , ad — be = 1 



(C23) 



a + 6r ' 

the doubly periodic elliptic functions with different periods are expressible through each 
other. 

Any transformation ( 10231) can be represented as a product of powers of two generators 



S := 



1 
-1 



T := 



1 
1 1 



(C24) 



The matrix S, generates the so-called Jacobi imaginary transformation 



T ^ t' := —1/t 



under which the Jacobian functions vary as 
sn[tz,fi)=t — -, cn[tz,fi) 



cn{z, jj) ' 



1 A f >\ dn(z,/i) 

— -, dn{iz, /i ) = — 

cn[z,jj) cn[z,jj) 



(C25) 



where /x' := yl — /i^, the complementary modulus. The action of the generator T which 
represents the shift transformation 

r ^ r' := 1 + r, 



We assume that Re(T) > . 
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results in the following relations between the basic Jacobian functions 

^"(^"V^^^d^^' ^"^^"V^^d^^' ^"^^"V^^d^^- ^^26) 

Finally the change 

t^t' :=t/(1 + t), 

which can be represented as, TST, is of special interest to us since it relates the functions 
with the inverse moduli /i and := 

sn{ij,z , n) = n sia{z , fj.) , cn{f^z , ii) = dn{z , fi) , dii{nz , ii) = cn{z , fj.) . (C27) 
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